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ABSTRACT 

We have performed hydrodynamical simulations from cosmological initial conditions 
using the AMR code RAMSES to study atomic cooling haloes (ACHs) at z = 10 
with masses in the range 5 x 1O 7 M <J M <;2 x 1O 9 M0. We assume the gas has pri- 
mordial composition and H 2 -cooling and prior star-formation in the haloes have been 
suppressed. We present a comprehensive analysis of the gas and DM properties of 19 
haloes at a spatial resolution of ~ 10 (proper) pc, selected from simulations with a to- 
tal volume of ~ 2000 (comoving) Mpc 3 . This is the largest statistical hydro-simulation 
study of ACHs at z > 10 to date. We examine the morphology, angular momentum, 
thermodynamical state, and turbulent properties of these haloes, in order to assess the 
prevalence of disks and massive overdensities that may lead to the formation of super- 
massive black holes (SMBHs). We find no correlation between either the magnitude 
or the direction of the angular momentum of the gas and its parent DM halo. Only 
3 of the haloes form rotationally supported cores. Two of the most massive haloes, 
however, form massive, compact over-dense blobs, which migrate to the outer region 
of the halo. These blobs have an accretion rate ~ O.5M yr _1 (at a distance of 100 
pc from their center), and are possible sites of SMBH formation. Our results suggest 
that the degree of rotational support and the fate of the gas in a halo is determined 
by its large-scale environment and merger history. In particular, the two haloes that 
form over-dense blobs are located at knots of the cosmic web, cooled their gas early on 
(z > 17), and experienced many mergers. The gas in these haloes is thus lumpy and 
highly turbulent, with Mach numbers M. > 5. In contrast, the haloes forming rota- 
tionally supported cores are relatively more isolated, located midway along filaments 
of the cosmic web, cooled their gas more recently, and underwent fewer mergers. As a 
result, the gas in these haloes is less lumpy and less turbulent (Mach numbers M<4), 
and could retain most of its angular momentum. The remaining 14 haloes have a di- 
verse range of intermediate properties. If verified in a larger sample of haloes and 
with additional physics to account for metals and star-formation, our results will have 
implications for observations of the highest-redshift galaxies and quasars with JWST. 

Key words: galaxies: formation — large-scale structure of the universe — stars: 
formation — turbulence. 



1 INTRODUCTION 

According to the paradigm of hierarchical structure forma- 
tion, the dissipationless dark matter (DM) particles build up 
structures in the Universe from small building blocks to cre- 
ate increasing more massive DM haloes. Due to its ability to 
cool radiatively, the baryonic matter can loose pressure sup- 
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port and condense inside the DM haloes, forming the first 
luminous objects: stars, galaxies, and black holes (BHs). 

The formation of the first stars has been extensively 



studied in the last decade by a number of authors (e.g. Abel 
|etaL] POOP) ; | Yoshida et aT] p006t ; [Turk et al. | ( [2012] ) ) . These 
studies have shown that the first stars (pop III stars) were 
formed due to molecular hydrogen cooling inside DM haloes 
of ~ 10 6 Mq at redshift as high as z ~ 30. Numerical simu- 
lations of the first stars were extended to the DM haloes in 
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the mass range ~ (10 7 - 1O 8 )M DM by e.g. 



Wise & Abel 



i 2007[ ) , |Greif et al.| ( |2008[ ) , and |Prieto et al.| ( |2012[ ) in order 
to study the physical processes involved in the formation of 
the first galaxies. Unfortunately, individual pop III stars (or 
even clusters of hundreds of these stars) are too faint to be 
detected even with current (and future) observations. 

The possibility that the James Web Space Telescope 
(JWST) can observe the first galaxies, located in haloes 
J>10 times more massive than those hosting the first PopIII 
stars, is an exciting motivation to study the so-called atomic 
cooling haloes (i. e. DM haloes with a virial temperature 
> 10 4 K; hereafter ACHs - see e.g. Bromm & Yoshida 
(20111; Ferrara (20101 for a recent reviews). These types ot 



haloes have been studied recently by |Romano-Diaz et al.| 
( |2011[ ) and |Pawlik et~aL] ( |2011[ ) in the context of galac- 
tic disk formation with smoothed particle hydrodynamic 



(SPH) simulations. Romano-Diaz et al. (20111 presented a 



low-resolution statistical study of haloes with masses above 
~ 10 8 Mq inside a highly biased over-dense region at z = 10, 
showing that every halo with mass above ~ 1O 9 M0 devel- 



ops a gas disk structure. Pawlik et al. (20111 have stud- 
ied the formation of a galactic disk inside a sing le 1O 9 M 
halo using a high-resolution SPH simulation at the same 
redshift. Wise et al. (20121 similarly studied the formation 



of two proto-galaxies in high-resolution AMR simulations 
that also included star-formation and chemical and radia- 
tive feedback. Protogalactic disk formation inside ACHs has 
been studied using analytical (semi-phenomenological) ap- 



proaches as well. Under the simple assumption (Mo et al 



1998) that the disk specific angular momentum is a fixed 
fraction of the DM halo specific angular momentum, |Oh fc| 
Haiman ( 2002 1 found that most haloes are able to form a 



stable disks (the exceptions are haloes with very low spin 
parameter). Stringer & Benson (20071 extended the galac- 



tic disk study, taking into account the hierarchical nature 
of the structure formation. Under the assumption of con- 
servation of specific angular momentum during the collapse 
process they studied the effect of mergers on disk forma- 
tion, concluding that mergers have a significant impact on 
the prevalence of disks. 

Pristine ACHs have been proposed as possible sites for 
the formation of supermassive black holes (SMBHs). If the 
gas were able to remain at a temperature of ~ 10 4 K in these 
haloes, it has been proposed that fragmentation may be 



avoided, allowing the formation of a ~ 10 



M 



q supermas- 



Koushiappas et al.| 


2004| |Volonteri & Rees| 2005| Lodato 


fc Natarajan||2006 


Spaans & Silk 2006| Begelman et al. 



2006j). These SMBHs could then grow by another factor of 



10 



3=4 



by redshift 



6 and explain the puzzling presence 



of M J>10 9 Mq BHs by those redshifts (see, e.g., the recent 
Haiman|20f2] ). A few groups have now performed 



review m 



simulations of gas collapse i n such haloes (|Wise et al. 2008 
Regan fc Haehnelt|2009a|b||Latif et al.|2011[ ), with the sim- 
plifying assumption that H2 cooling can be avoided, which 
may be feasible in case the haloes are illuminated by a strong 



Lyman Werner background (Omukai 2001 Bromm & Loeb 



2003 Shang et al.|2010[ ). These works have indeed found no 
evidence of fragmentation. 

In this paper we continue the study of ACHs using AMR 
simulations. Our work is motivated by the possible forma- 
tion of galactic disks and SMBHs inside these haloes, as well 



as by the fact that these haloes may host the first galaxies 
that will be directly detectable by JWST. Our study is the 
largest in the literature based on hydrodynamical simula- 
tions, with a total of 19 haloes analysed at a spatial resolu- 
tion of Ax ~ 10 pc. We present, for the first time, a compre- 
hensive study of gas and DM physical properties, including 
their spatial distribution, angular momenta, thermodynami- 
cal state, and turbulent properties. We go beyond commonly 
used diagnostic tools, such as spherically averaged profiles, 
and compute density and velocity power spectra, and den- 
sity, velocity and temperature probability distribution func- 
tions. Using our sample of haloes, we attempt to identify the 
physical mechanisms that determine the angular momentum 
of the gas in the inner halo, and therefore whether a rota- 
tionally supported inner disk or a low-spin, compact dense 
cusp (and ultimately a SMBH) forms. 

Although we find essentially no correlation between the 
magnitude or direction of the angular momenta of the gas 
and the DM haloes, there are clear differences between the 
haloes hosting the two types of inner gas: haloes with low- 
spin central cusps form in highly biased regions at super- 
knots of the cosmic web, and have cooled their gas early 
on and undergone many mergers. This results in a lumpy 
and highly turbulent medium that could produce efficient 
transfer of angular momentum. On the other hand, the ro- 
tationally supported core is found in haloes that are located 
midway along cosmic filaments, cooled relatively more re- 
cently, and have experienced fewer mergers. The gas in these 
haloes is less lumpy and less turbulent and could retain most 
of its angular momentum. By a detailed study of the turbu- 
lent properties of the intra-halo medium (hereafter IHM) we 
have found that turbulence plays a major role at dissipat- 
ing (or preserving) angular momentum via fragmentation of 
the IHM. Turbulence in turn is fed by the merging process. 
We will show that compact over-dense blobs (COBs) form 
in highly biased regions at the centers of super-knots, while 
rotationally supported cores (RSCs) form in the middle of 
filaments in nearly un-biased regions. Our numerical experi- 
ments suggest that neither rotationally supported disks nor 
SMBHs are a common occurrence in ACHs, and should be 
followed up by studying a larger sample of haloes, with more 
realistic gas physics and more detailed analyses of turbulent 
angular momentum transport. 

The paper is organized as follows: In section ^2] we de- 
scribe the details of our numerical simulations and our halo 
sample. In Sj3] we explain our analysis and show the physical 
quantities extracted from the simulations. In section ij4] we 
discuss our results and finally in i|5]we present our conclu- 
sions. 



2 NUMERICAL SIMULATION DETAILS 

We used the AMR code RAMSES ( |Teyssier||200"2"l> with a 
modified non-equilibrium cooling module ( |Prieto, Infante, fc 
Jimenez|2008 1 in order to study the main physical properties 
of both DM and gas components of 19 ACHs in the mass 
5 x 10 7 M Q to ~ 2 x 10 9 M Q . We assumed the gas 



range 

has primordial composition, and switched off all molecular 
cooling processes, allowing the gas to cool down only due to 
atomic hydrogen and helium lines. 

We worked in a concordance ACDM cosmological 
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model: h = 0.719, 
as — 0.796 and n s 



D.A = 0. 742, n m = 0.258, fl b = 0.0441, 
0.963 ( [Komatsu et~ah|2011[ ). The dy- 



namical ini tial conditions were taken from mpgrafic ( Prvmet 
et al. 1 |2008[ ) and the initial abundances of H, H + , He, He 



Galli & Palla 



(1998) at 



and He ++ , D, D + and e~ from 
z — 100 (the initial redshift of the simulations). 

First, we performed a number of DM-only simulations 
with 256 3 particles inside (3 Mpc) 3 , (5 Mpc) 3 , (5.5 Mpc) 3 
and (8 Mpc) 3 boxes in order to look for ACHs, i.e., haloes 
with virial temperatures T v i r > 10 4 K, which corresponds to 
a Af v ir > 10 8 M Q at z = 10. The relation T vir (M vir ) was taken 
from Bryan & Norman (1998) and M v i r refers to the mass 
for a DM overdensity with a density contrast 8 = 200. In 
this study, we focus on z = 10, which corresponds to the 
epoch when ACHs are "2-er" objects and begin to have a 
significant abundance (a few per Mpc 3 ). Using the HOP al- 
gorithm ( jEisenstein fc Hut|1998[ ) we identified DM haloes in 
the DM simulations outputs at this redshift. Our goal was to 
study the simplest, relatively isolated first protogalaxy, and 
we therefore avoided haloes that are undergoing an obvious 
merger at z — 10. We have found that mergers are very com- 
mon at high redshift, and we have rejected many simulations 
(containing ~ 100 pre-selected haloes) in which no such iso- 
lated DM halo was found, before arriving at a sample of 19 
suitable, non-merging DM haloes. We then performed a new 
set of gas+DM simulations for these 19 haloes, in 19 different 
runs. Each of these re-simulation runs were performed with 
512 3 DM particles inside the corresponding original simu- 
lation volume, from z = 100 to z = 10 in order to follow 
the main process of gas collapse into the ACHs at higher 
resolution. In Table [l] we summarize the properties of our 
19 haloes, and the simulations in which they were identified. 
The virial quantities are computed from the halo mass tak- 



ing the expressions of Bryan & Norman ( 1998 1. Each of our 
simulations used ~ 60, 000 CPU hours (for a total of ~ 10 6 
total CPU hours for the whole simulation suite) on a cluster 
with 500 cores and 2 TB of distributed memory. 

Throughout this study, we adopted the following refine- 
ment strategy: if there were more than 4 DM particles inside 
a cell, then the cell is refined; if the gas density inside a cell is 
more than 4 times the average density inside the whole box, 
then the cell is refined, and if the linear cell size is larger 
than 0.25 its corresponding Jeans length, then the cell is 
refined ( |Truelove et al.|1997[ ). The maximum level of refine- 
ment was l m ax = 15, corresponding to a proper resolution 
at z = 10 of 8, 14, 15, and 22pc in the (3 Mpc) 3 , (5 Mpc) 3 , 
(5.5 Mpc) 3 and the (8 Mpc) 3 box, respectively. Our resolu- 
tion corresponds to ~ 1% of the virial radius and it is not 
adequate to resolve structures in the inner core of the halo. 
We therefore set a gas temperature floor, Tfl oor , in order to 
avoid numerical fragmentation. This temperature floor was 
based on the Truelove criterion that the Jeans length must 
be resolved by 4 cells, 



Aj = 4Ax, 



(1) 



where Aj is the Jeans length, Axi is the proper cell size at 
a given level of refinement I. 

In their recent study of primordial star formation study 
with magnetic field and non-equilibrium H2 chemistry, |Turk| 
et al. (20121 show that a much higher number (64) cells per 
Jeans length is required to achieve convergence in the ultra- 
dense center of the halo. Since here we do not include H2 
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Figure 1. Baryon mass fraction inside different radii as a function 
of the halo mass. The horizontal line shows the average cosmic 
value of ni,/f! m =0.171 (using the parameters in the text). Circles, 
triangles, and squares show the mass fractions inside 0.5, 1, and 
2i? v ; r , respectively. 



chemistry or magnetic fields, and we are not interested in 
studying the gas collapse process at sub pc scales, the stan- 
dard Truelove criterion should be adequate for our purposes. 
We then obtain the gas temperature floor at a given level of 
refinement as 



tfloor — 



64-7T Gfinipp L 2 



15 k E 



2 2( ' 



(2) 



where G is the Newton's gravitational constant, /1 is the 
mean molecular weight, m p is the proton mass, p is the gas 
mass density, ks is the Boltzmann constant and L is the 
proper simulation box size. Assuming p — 1.22 (appropriate 
for the nearly neutral T <C10 4 K primordial gas in our case) 
and I — lmax, Tfloor at z — 10 can be written as 



tfloor 



60K 



V 1 cm -3 



L 



5.5 Mpc 



(3) 



This temperature floor is much lower than the atomic cool- 
ing temperature threshold ~ 10 4 K, so it does not affect the 
thermal evolution of the gas unless its number density ex- 
ceeds > 100 cm -3 . In most of our haloes, the gas does not 
reach such high densities and the temperature floor is not ac- 
tivated. However, there are two exceptions, where the num- 
ber density in the core reaches up to ~ 10 3 cm -3 . In these 
cases the gas temperature in the central over-dense regions 
does not have a physical meaning (nevertheless, as we will 
see below, our results are sufficient to argue that conditions 
in these cores may be most favourable for forming SMBHs). 



3 RESULTS 

In this section, we present an exhaustive list of physical 
quantities extracted from our 19 simulated haloes. These 
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Table 1. Properties of the 19 atomic cooling haloes at z = 10 analyzed in this paper, and the simulations in which they were identified. 
The first 4 columns indicate the id number, mass, virial temperature, and virial radius of each halo, ranked by its total mass. The 5 th 
column shows the type of gaseous object forming in the halo (rotationally supported core or compact over-dense blob). The 6-8 th columns 
show, respectively, the comoving box size, DM particle mass, and spatial resolution for the gas (in proper units, at the highest level of 
refinement at z = 10). 
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include traditional global measures (such as gas mass frac- 
tions, as well as the geometrical shapes and spin parameters 
for both the DM and the gas), radial profiles (density, tem- 
perature, and angular momentum), three-dimensional mor- 
phologies, as well as phase diagrams for the gas. In addition, 
we present several less commonly used diagnostics, such as 
probability distribution functions (gas density, temperature 
and velocity), and a Fourier analysis (including power spec- 
tra of gas density and velocity), which will help us to un- 
derstand the physical conditions in ACHs at their collapse 
at z = 10. 



3.1 Gas fractions and phase diagrams 

We computed the baryon mass fraction ft, = M gas /(MDM + 
M gas ) inside each of the 19 haloes at three different radii: 0.5, 
1, and 2i? v i r . Here the virial radius 7? v j r is measured from the 
center of mass of the DM halo. Figure [T] shows these mass 
fractions as a function of the halo mass. The average baryon 
mass fractions are (/i,,o.5, /&,i.o> /f>,2.o) = (0.066,0.097,0.16), 
which correspond to 39%, 57% and 92% of the cosmic av- 
erage value (fb — 0.171). This is lower than the observed 



average value found by Romano-Diaz et al. (20111 at the 
same redshift and for the same range of masses using SPH 
simulations. Our ft, value is also lower than values typically 
found in galaxy cluster simulations at low rcdshifts (e.g., 
Gottlober & Yepes|(|2007|),|Ettori et al.|(|2006|) andlKravtsov 



et al. (20051) 



due to the fact that the gas temperature is not too far be- 
low the virial temperature, so that gas pressure can affect 
the dynamics and delay the collapse of the gas. Interest- 
ingly, however, the mass fractions lack any obvious trend 
with halo mass (or with T v i r ), as one would expect. The lack 
of this trend suggests that turbulence, generated at the time 
of virialization (Wise fc Ab"eT|2007[ ), can provide additional 
pressure support, and delay gas collapse even in more mas- 
sive haloes, where thermal pressure is less important. The 
hypothesis that gas has not yet fully collapsed in the haloes 
will be supported by our results below on spin misalignment 
angles and the gas spin parameters. 



Harford et al. ( 2008[ l show that the baryon mass fraction 



The low baryonic mass fractions we find might be partly 



increases inside smaller radii and it reaches a value near the 
cosmic average at R v i r . The ft, values in Figure [l] decrease as 
the average volume (and radius) decreases. In other words, 
there is no evidence of an overall concentration baryons near 
the DM halo center of mass in our simulations, suggesting 
that thermal and turbulent pressure can delay gas collapse 
even inside 0.5i? v ir- 

Figure [T] also shows a large scatter in ft, with halo-to- 
halo variations by an order of magnitude. Despite the low 
average fb at i? v ir, we find three haloes with /{, near the 
universal value even at 0.5i? v ir- As we will show below, these 
haloes are the ones with high central gas spin parameters 
and low misalignment angles between gas and DM spin. 

In Figure [2] we show the density-temperature n — T 
phase diagrams for gas inside iivir- The figure shows that 
at the highest densities, most of the gas is heated by shocks 
and by compression above the ~ 10 4 K limit, but radiative 
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Figure 2. Phase diagrams for the gas inside the virial radius in the 19 simulated haloes. 



cooling keeps the gas temperature near this value. At lower 
densities, a range of temperatures is seen, with a significant 
amount of cooler gas that has not yet been heated and some 
gas still significantly above 10 4 K. Also seen in this figure 
are a few haloes in which the gas density is prevented from 
reaching high values (n <Tcm -3 ). As we will see below, some 
of these haloes have a high degree of rotational support near 
their center (including three haloes, #1, #7, and #12, with 
full rotational support). 



3.2 Shape of the DM halo and its gas content 

In order to quantify the morphology of both the gas and the 
DM components, we computed ratios of the three semi-axis 
c < b < a of the ellipsoids, characterized by the eigenvalues 
of the gas/DM moment of inertia tensor inside the virial 
radius J? v ir- The moment of inertia tensor was computed 
from its standard definition 



the 



m a (r a 5ij 



(4) 



where i and j lab el Cartesian coordinates x, y and z, 
r = y^rl + Ty + r\ is the distance from the DM (gas) cen- 
ter of mass to the a th DM particle (gas cell) and Sij is the 
Kronecker delta function. 



semi-axis follow the relations 



(5) 



'Ai b = Xi 
A 3 ' a \ \ 2 ' 

Figure |3|a) shows the semi-minor to semi-major axis 
ratio c/a as a function of the halo mass. It is clear that the 
DM component (squares) tends to be more spherical than 
the gas component (circles). The average value is c/a = 0.62 
for the DM and c/a — 0.47 for the gas. Figure |3^b) shows 
the fraction of haloes with the gas c/a ratio below 0.5. For 
masses above ~ 6 x 1O 8 M0 all haloes show c/a < 0.5. This 



plot is similar to Figure 3(c) in Romano-Diaz et al. (20111, 



from which they conclude that the disk fraction is 100% for 
more massive haloes. However, as we will see below, despite 
of the low c/a, it is far from being rotationally supported. 

In general, the set of criteria that defines a conventional, 
rotationally supported disk should include: (i) low c/a ratio 
together with b/a ~ 1; (ii) alignment of the gas angular 
momentum vector with shortest principal axis c and (iii) a 
gas spin parameter A w 1. Figure [4] shows the 19 haloes 
in the c/a — b/a plane. As this figure illustrates, the DM 
component (squares) tends to populate the high c/a ratio 
region in the plane, while the gas component (circles) tends 
to reach much lower values. This would be expected if the 



6 



Prieto, Jimenez & Haiman 



(b) 



□ □ 

o 

□ 

G 




0.95 - 



0.9 



0.85 



0.8 



0.75 



0.7 



0.2 



#12 



□ □ 
□ 



DM □ 

3as 



oo 



m 



o □ 
o q 



#7 □ 



□ □ □ 



0.5 

c/a 



Figure 3. (a): The ratio between the semi-minor and semi-major 
axis (c/a) for both gas (circles) and DM (squares), as a function 
of halo mass. The gas component tends to have lower c/a than the 
DM component, (b): The fraction of haloes with gas component 
c/a ratio below 0.5. Despite of the fact that the fraction of haloes 
with c/a < 0.5 is 100% above ~ 10 9 Af Q , we find no rotationally 
supported galactic discs above this mass limit. 



Figure 4. The gas (circles) and DM (squares) principal axis ra- 
tios in the c/a — b/a plane. The DM occupies the high c/a, b/a 
region of this plane, while the gas component tends to occupy 
lower values: the DM is more spherical than the flattened gas 
component. Both components typically have a significant triaxi- 
ality. 



gas is able to collapse in the direction perpendicular to the 
angular momentum vector to assume a more flattened shape. 
We computed the misalignment angle between the angular 
momentum and the shortest semi-axis "a" (inside 7? v ir) and 
found that there is an offset of between 30-90 degrees in 
the whole sample. Overall, the shapes of the DM and the 
gas are both typically quite triaxial. The shapes of the gas 
and DM nearly remain closer to prolate rather than oblate 
[c/b < b/a), with average b/a ratios of b/a = 0.84 for the 
gas and b/a = 0.87 for the DM. This feature has been seen, 
e.g. in the DM-only simulation of Vera-Ciro et al. (20111 



where the DM matter haloes develop a prolate shape at the 
beginning and evolve to an oblate shape at later times. 



In order to study the global rotational features of our 
systems, we computed the spin parameter A' defined by Bul- 
lock et al.|(|2001a|): 



a: 



j, 



gas J DM 



gas/ DM 



s/DM 



(7) 



The definition of A' in eq. [7] corresponds to \/\f2 times 
the ratio between the angular velocity ui and the angular 
velocity ui a of a rotationally supported system with the same 
mass and radius (i.e. M gas / DM , the gas/DM mass inside the 
virial radius). J gas /DM 



J. 



gas/DM 



E 



\Jgas/DM\ with 



WaT ' a ^gas / 'DM X V a ^ gaa / DM 



(8) 



3.3 DM halo vs. gas angular momentum 



Doroshkevich (19701; White (19841; Heavens & Peacock 



( 1988 1 showed that the initial angular momentum of a sys- 



tem in a hierarchical Universe is proportional to both the 
tidal tensor Tn and to the moment of inertia tensor Lj 



Ji oc eijkTjilu 



(6) 



The evolution of the angular momentum of DM haloes in a 
cosmological context has subsequently been revisited by. e.g. 
|Vitvitska et al.| ( |2002[ ) , to study the importance of mergers. 
They found that the angular momentum typically has in- 
creases sharply in major mergers, and declines steadily dur- 
ing periods of gradual accretion of small satellites. The an- 
gular momentum evolution has moved from a picture where 
the tidal forces due to the neighboring structures are solely 
responsible for the halo spin to a picture in which the merger 
history of the halo has a significant impact. 



is the gas/DM angular momentum. This quantity is com- 
puted with respect to the gas/DM center of mass and the 
velocities v gaa / DM are corrected by the gas/DM center of 
mass velocity; m a is the gas/DM mass at position f a . All 
these quantities are computed inside the virial radius _R v i r 
of the halo. 



It has been shown by a number of works (e.g. Bullock et 
al. (2001a I and references therein), that the spin parameter 
A follows a log-normal distribution: 



hi(A')]dmA' 



MA'/A') 



din A 



(9) 



where A' is defined by In A' = lnAi)/iV (the geometric 

average of the spin parameter for N 3> 1 haloes) and a is its 
standard deviation. |Bullock et al.| ( |2001a[ ) have shown these 
parameters have values A' « 0.04 and a « 0.5. 
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In order to study the correlation between the DM and 
gas spin orientations, we also computed the misalignment 
angle 9 between their spin directions, 



cos(9) 



(10) 



Figure [5^a) shows the spin parameters for both gas 
(circles) and DM (squares) as a function of halo mass. 
The average values over the 19 haloes are nearly identi- 
cal, Arj M = 0.046 and A gas = 0.047. The average DM spin 
value is in good agreement with that found in a previ- 
ous study using N-body simulations and focusing on high- 
redshift haloes (Davis & Natarajan 2009). Overall, the result 
Ajj M ~ Ag as = 0.047 is also consistent with the usual simple 
picture in which the DM and gas initially have similar spin 
parameters at the virial scale. However, we do find an aver- 
age ratio of (A(, aa / Af) M ) = 1.65. This value is h igher than the 
values reported by |Gottlober fc Yepes| ( |2007| > for > 10 14 M Q 
haloes at low redshirt, and also higher than the simple ratio 
(Agas)/(ArjM) — 1.02. This is because the gas spin param- 
eter tends to be noticeably smaller (by a factor of ~two) 
than the DM spin parameter at masses Mh<lO 8 M0. At 
higher masses, the gas spin parameters increase and exceed 
the DM spin in some cases. This tendency is illustrated more 
clearly in Figure |5|b), where we show the ratio between the 
gas and DM spin parameter as a function of the halo mass. 
Figure [5^c) shows Ag as v/s Ad M . There is no a clear trend 
among these two quantities. This fact will be supported be- 
low by the measurement of the misalignment angle among 



gas and DM angular momentum. |Kimm et al. (20111 have 



recently studied the relative angular momenta of gas and 
DM haloes, and also found significant differences. In partic- 
ular, they found that the discrepancy between the specific 
angular momentum of the gas and the DM in the region 
0.1 < r/i? v ir < 1.0 is due to the cooling process governing 
the gas dynamics. The gas is accreted through cold filaments 
into the inner regions of the halo, where it condenses; how- 
ever, the accreted DM crosses the central region and redis- 
tributes the angular momentum in the outer regions of the 
halo. 

Figure [6] shows the angle between the DM and the gas 
angular momentum vectors. The DM and gas spins are mis- 
aligned for all 19 haloes we analyzed, by an average angle 
(6) = 91°. This angle is almost 5 times the average misalign- 
ment angle of around 20° found by Sharma & Steinmctz 
et al. (20051 at z = 0, using the SPH technique. Figure [6 
shows no trend of 9 with halo mass, and clearly shows that 
despite their similar magnitudes, the gas and DM spins must 
have different physical origins. Interestingly, we find that the 
three most aligned systems have high gas spin parameters 
and high baryonic mass fractions. This suggest that the gas 
in these haloes has collapsed enough to be aligned. 

Figure [7] shows the correlations between the baryonic 
mass fraction fb{< 0.5it v ir), and the gas spin parameter 
A gas within both < lR v i r (squares) and < 0.5i?vir (circles). 
There is a clear correlation between A gas (< 0.57? v ir) and 
fb(< 0.5-Rvir), showing that gas inside haloes with more col- 
lapsed gas reaches a higher spin. In particular, the three 
haloes with the lowest misalignment angles present high fb 
and Agas • This is consistent with the simple picture in which 
the DM and gas have similar angular momenta at larger 
radii, but the gas can continue to collapse along their mu- 



st 1 
o 



<< 1(T 



5 

Q . 
h 10° 



<< 10-' 



□ 





DM □ 

3as 



+ 
+ 



+ 

h + 
+ + + 



2 

o 

k 



+ + * + ++ 

+ + + + +± + 

+ + * + 



(C) 



gas 



Figure 5. (a): The "global" spin parameters, measured within 
-Rvin for both gas (circles) and DM (squares) as a function of 
halo mass. The gas component has lower spin than the DM haloes, 
with a tendency to match at higher masses, (b): the ratio between 
the gas and DM spin parameters as a function of halo mass. The 
ratio between the spin parameters tends to increase at higher 
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Figure 6. Misalignment angle between the DM and gas spin 
vectors as a function of halo mass. There is no clear trend for the 
misalignment with halo mass. This indicates that gas and DM 
angular momenta have different physical origins, likely influenced 
by the local environment rather than determined solely by tidal 
torques (see text for more discussion). 
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Figure 7. Spin parameters A gas at < R v i r (squares) and < 
0.5-R v i r (circles) as a function of the baryonic mass fraction 
within 0.5R v i r . There is a very strong correlation between the 
spin parameter and within 0.5-Rvin showing that the gas inside 
haloes with more collapsed gas eventually reaches higher spins. 



tual angular momentum axis, and is spun up by a larger 
factor - eventually this could lead to disk formation. 

3.4 Radial profiles 

We computed the radial profiles of different quantities for 
both the gas and the DM component of each halo. This in- 
cludes the mass-weighted specific angular momentum and 
the mass density, separately for the gas and the DM. For 
the gas component, we additionally computed the mass- 
weighted temperature. All quantities were computed as 
spherical averages centered at the component (gas/DM) cen- 
ter of mass. Note that this, in general, does not coincide with 
the point with the highest density. 

More specifically, the mass-weighted average value q of 
a physical quantity q was computed as 

q{n)m{n) 



(f) 



Mr 



(11) 



where r; is the position of each DM particle or hydro cell in- 
side an spherical element of volume given by AV = 4irf 2 Ar. 
The averaging bin Ar takes the value Ar = 8Ax for gas and 
Ar = 16Ax for DM, with Aa; the minimum resolved scale 
of the simulation; m(ri) is the DM or baryonic mass at the 
position Vi. Mr is the total DM or baryonic mass inside the 
averaging volume and f = r + Ar/2. 

Density profiles. Figure|8]shows both the gas (solid line) 
and DM (long-dashed line) density profiles. As mentioned 
above, the DM and gas component are centered at their re- 
spective center of mass. This is the reason why the highest 
densities are not located at r — 0. For reference, the short- 
dashed line shows the isothermal radial profile p oc r~ 2 and 
the dotted line show a p oc r -3 , the NFW density profile at 
large radius ( jNavarro, Frenk'fc White 1997}. These power- 
law profiles are not good descriptions of the majority of the 



profiles (especially for the baryons) but are good approxi- 
mations over limited ranges of radii (especially for the DM 
profiles in the outer halo near i? v ir). 

Temperature profiles. The temperature profiles are 
shown in Figure [9] As expected from ACHs, the gas tem- 
perature is almost constant around ~ 10 4 K inside the virial 
radius. This temperature is not affected by the numerical 
temperature floor in eq. (Jsj, because at the densities we re- 
solve inside these haloes, Tfl oor is much lower than 10 4 K. 
The two exceptions are the sharp peaks in the temperature 
profiles of haloes #16 and #18 near R v i T (in the bottom 
row). These correspond to strong overdensities (as shown in 
Fig.[8| which formed at z ~ 11. Despite our poor spatial res- 
olution, as we will argue below, these runaway overdensities 
are candidate sites for SMBH formation. These objects have 
masses of M b i ob Ri2x 10 7 M Q and M blob w 3 x 1O 7 M for 
haloes #16 and #18 (with halo masses of 9.22 x 1O 8 M and 
1.53 x 1O 9 M , respectively). Interestingly, in both cases, the 
strong gas overdensities are located near the virial radius, 
rather than near the halo center. It is tempting to conclude 
that we are seeing a merger, and to associate the overdensi- 
ties with the cores of the merging partner haloes. However, 
this is not the case: the overdensities are seen only in the gas 
and not in the DM. In other words, the dense gas "blobs" are 
not at the center of another DM halo. We have examined the 
evolution of the gas distribution over time, and have found 
that the overdensities are produced in large-scale, highly su- 
personic shocks, and are originally located near the center 
of the halo. The gas flows producing the blobs are not per- 
fectly radial toward the center of the halo, however, and the 
interacting gas shows a velocity relative to the halo's center 
of mass. The over-dense blob itself, produced by the strong 
supersonic shocks, has a significant residual velocity relative 
the halo's center of mass. This residual velocity exceeds the 
escape velocity from the central region of the halo, allowing 
the blob to travel to the outer region of the halo within a 
--^dynamical time. A similar phenomenon has been seen in 



numerical simulations of star- forming regions ( Padoan et al 
[19971 



To examine these blobs further, in Figure [To] we show 
their radial profiles at two different redshifts, re-centered 
at the center-of-mass of each blob. Row (a) shows the gas 
(number) density profiles, and row (b) shows the tempera- 
ture profiles. As mentioned above, inside r ~ 100 pc, the 
gas temperature increases unphysically, because the activa- 
tion of the temperature floor (we emphasize that our results 
are conservative, in the sense that without the temperature 
floor, the blobs should have collapsed further and reached 
even higher densities). Row (c) shows the enclosed gas mass 
and row (d) shows the specific angular momentum, divided 
by the Keplerian value at each radius. This ratio is always 
below unity, but increases toward smaller radius, reaching 
a maximum value of L/Lk ep ~ 0.5 at the innermost re- 
gions. This value is similar to the values shown by |Regan fc| 
Haehnelt (2009a I in their Figure 7 at few x 10 pc. As in that 



paper, these blobs can be similarly argued to be sites where 
compact disks, and eventually SMBHs, may form. Due to 
the lack of our spatial resolution, however, we can not say 
more about the inner structure of these blobs. 

We next compute the mass accretion rates (at z — 10) 
onto the dense blobs as follows: 
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Figure 8. Radial mass density profiles for gas (solid line) and DM (long-dashed line). For reference, the short-dashed lines and the 
dotted lines show power-law slopes corresponding to an isothermal profile (p <x r~ 2 ) and the NFW profile at large radius (p <x r~ 3 ). The 
profiles of the most massive haloes reveal the presence of over-dense blobs at ~ R v ir> which we argue are candidate sites for the formation 
of SMBHs. These blobs result from large-scale interacting shocks near the halo center. Due to the asymmetry of the interaction, the 
blobs retain a non-zero velocity and are able to migrate to the halo outer region (see text). 



M b , 



= 4irr 2 ( piVi 



(12) 



where r*j is the gas grid position measured from the blob's 
center of mass, r% its modulus, and f» is the unit vector in 
that direction; f is the radius averaged in a finite shell, pi 
is the density at fi and Vi is the gas velocity with respect 
to the blob's center-of-mass velocity at fi. We have found 
Mttiob ~ 0.5 MQyr -1 at a distance r w 100 pc. This value 
should, of course, be taken only as an upper limit for the 
rate that would reach any putative SMBH forming near the 
center of the COB. This rate is roughly in agreement with 
the rates ~ M©yr -1 typically found inside atomic cooling 



haloes ( [O'Shea fc Norman|[2007l |Shang et aLl[2010| . This 
latter value is of order the mass accretion rate ~ ci/G oc 
T s/2 /G expected in a self-gravitating gas at temperature 
~ 10 4 K. 

Angular momentum profile. Figure [TT] shows the pro- 
files of the specific angular momentum (SAM) for both gas 
(in solid line) and DM (in long-dashed line). Overall, the 
DM and gas SAM trace each other well at small radii in 
most of the haloes. This agreement break down outside 
~ -Rvir, where the DM SAM exceeds the gas SAM, and in- 
creases more rapidly with radius, as well. The short-dashed 
line shows the Keplerian SAM profiles, computed for the 
density profiles found in the simulations (Fig. |8p . At radii 
r > 0.3-Rvir, the Keplerian curves are higher than the real gas 



SAM, showing that the gas is not rotationally supported. 
At smaller radii there are three haloes (#1, #7, and #12, 
with Mum = 5.354 x 1O 7 M , Mum = 3.987 x W 8 M Q and 
Mum = 5.718 x 10 8 M Q ) in which the gas SAM reaches the 
Keplerian value, suggesting the existence of rotationally sup- 
ported disks in the inner regions of these three haloes. How- 
ever, inspection by eye does not reveal geometrical disks for 
these three haloes. In the inner 0.1i? v ir, we find axis ratios 
of c/a ~ 0.6 and b/a « 0.8. In particular, the semi-minor to 
semi-major axis ratio c/a is close to the ratio c s /-u c i rc « 0.7 
between the sound speed and the circular velocity; consis- 
tent with thermally supported "fat disks" expected when the 
gas temperature is not significantly below the virial tem- 



perature (Oh fc Haiman 2002 Regan fc Haehnelt 2009a) 



Hereafter we refer to these objects generically as having ro- 
tationally supported cores (RSCs). Figure [l2"| shows the cor- 
responding gas and DM spin parameter profiles. As shown 
in Figure [5jb) , the gas spin parameter is higher than the 
DM spin parameter in some of the haloes at R vlr ; overall, 
there seems to be little correlation between the details of 
the two. For most of the haloes, Ag aa increases at smaller 
radii. As seen already in Figure [TT] there are three cases in 



which the gas spin parameter reaches Ag; 



1, this hap- 



pens at radii r<0.2i? v ir- The cores of these three haloes 
within this radius are rotationally supported. Interestingly, 
even in these three cases, there is a significant misalignment 
between the direction of the shortest (semi-minor) axis and 
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Figure 9. Mass-weighted gas temperature profiles. Due to atomic line cooling the gas temperature is almost flat at 10 4 K in the central 
regions. The peaks in temperature in two of the most massive haloes are associated with the overdensities in Fig. [8] and arc unphysically 
large (due to the activation of a numerical temperature floor). 




Figure 10. Radial profiles at the two overdensities at two different redshifts, including the number density (a), the temperature (b), the 
gas enclosed mass (c) and the specific angular momentum normalized to the Keplerian angular momentum (d). The first two columns 
correspond to the over density hosted by the A/dm = 9.220 X 10 8 Mq DM halo (#16) and the last two columns correspond to the 
A/pjM = 1-534 X 10 9 Mq DM halo (#18). The central temperatures are unphysically high because of the activation of the numerical 
temperature floor. Row (d) shows a maximum specific angular momentum ratio L/L^ep ~ 0.5. This value is similar to the value found 
by|Regan fc Haehnelt| (|2009a.[| at r fa fewxlO pc, supporting the notion that these blobs are promising sites of SMBHs formation. 
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Figure 11. Mass- weighted specific angular momentum profiles for the gas (solid curves) and DM (long-dashed curves). The gas and DM 
specific angular momenta trace each other and both decrease toward the center. The short-dashed line show the Keplerian specific angular 
momentum (SAM) for the measured density profile. At radii r>0.3H v i r the DM tends to have more specific angular momentum than 
the gas. At smaller radii, there are three cases with the gas SAM exceeding the Keplerian values, i.e. the gas in the core is rotationally 
supported. 



the angular momentum vector, by an angle 20° < <f> < 40° 
when measured inside 0.3-R v ir- 



3.5 Beyond spherically averaged profiles 

3. 5. 1 Power spectra of the intra-halo medium 

We next perform a Fourier analysis, in order to understand 
the properties of the turbulent environment through which 
the gas collapses into the DM haloes. The power spectrum 
of a physical quantity is the square of the Four ier transform 
averaged on spherical shells of radius k — k% + 
in wave number space. Following the above definition we 
computed the gas mass density p power spectrum (PS), in 
a box of side s ~ 4R v [ r around the halo's center of mass, as: 



P p = 47vk p k ■ pt , 

and the gas velocity PS 



K = -{-iirk 2 v k ■ vl), 



(13) 



(14) 



where Vk {pk) is the Fourier transform of the gas velocity 
(density), vl (pi) is its complex conjugate, and k = \k\ is 
the modulus of the wave number. For the velocity field, we 



have subtracted the motion of the center of mass, to avoid 
any spurious turbulence. The pre-factor of 1/2 in the veloc- 
ity power spectrum converts this quantity into the kinetic 
energy stored in wave numbers around k. With this factor 
included, we computed the fraction R k of the kinetic energy 
in solenoidal modes, 



Rk = 



|(47rfc 2 t3 a , fc ■ v^k) 
\(^k?v h ■ vl) ' 



(15) 



where the compressional (V x v c ,k = 0) and solenoidal 
(V ■ v s ,k = 0) components of the velocity field are given 
respectively by 



v c ,k = (vk ■ k)k/k 



and 



v s ,k — v k — (vk ■ k)k/k 



(16) 



(17) 



In turbulent fluid dynamics, there are two well-studied 
limits for the turbulent energy power spectrum: the k~ 2 ID 
Burgers spectrum and the fc _5//3 3D Kolmogorov spectrum 
(IKolmogorov | p4l) . The first one describes a system with 
velocity u and a fluid viscosity v governed by the equation 
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Figure 12. Gas and DM spin parameters corresponding to the SAMs in Figure [TT] 



du 
0t 



dx 



dx 2 



(18) 



For a high Reynolds number or v — > (applicable for this 



work), eq. 18 describes the evolution of a ID shock wave. 



Despite the one-dimensionality of the Burgers equation, it is 
often used as a reference for turbulent 3D fluids. A different 
3D treatment comes from the seminal work of |Kolmogorov| 
( |1941[ ). In this theory, the statistical properties of a turbu- 
lent fluid with viscosity v , mean energy dissipation rate per 
mass e and a minimum scale i associated to the energy dissi- 
pation (by viscosity) are fully determined by these three (or 
two if v = 0) quantities in the inertial scale range: L < I < I, 
with L the scale associated to the source of the turbulence 
(in our case L > i? v ir)- 

Figure [13] shows the density power spectra for the gas 
in the range of wave numbers corresponding to proper scales 



0.1 kpc <,2n/k ^6.3 kpc. Saichev & Woyczynski (19961 dis- 
cussed two limits of the scaling relation tor the density PS 
in the Burgers equation: P p oc k° and P p oc k~ 2 . The 
P p oc k° case corresponds to the limit of negligible pres- 
sure in which the fluid is dominated by strong shocks while 
P p oc k~ 2 arises in the presence of fluid pressure, i.e. a weakly 
shocked fluid. The gas density power spectra in our simula- 
tions do not show a well-defined single power-law, such as 



found in Gazol & Kim (20101. The PS is flatter on large 



scales, above ~ 1 kpc, whereas for most haloes, the power 
decreases steeply with k on smaller scales. The power-law 



slopes on large (sm all) scales are close to a as (a = —1.5). 
Kim & Ryu ( 2005 I have shown that the slope of P p flattens 



as the Mach number M increases. Following Gazol & Kim 
(20101, our power spectra can be interpreted as evidence for 



a turbulent cascade from highly supersonic turbulence (with 
large M) on the virial (shock) scales of ~kpc, to low-AI tur- 
bulence on small scales (higher k). 

Probably the most significant feature in these power 
spectra is that for two of the most massive haloes, the PS in- 
creases at large k. These behavior is characteristic of systems 
with self-gravity as shown by |Collins et al.| ( |2012[ ), where the 
densest regions have enough mass to collapse gravitationally 
and form gas clumps increasing the power on small scales. 
As mentioned by |Collins et al.| ( |2012[ ) , this feature in the 
density PS is indicative of a highly lumpy environment. In 
our case, this occurs for the two haloes which developed 
over-dense blobs, and clearly distinguishes these two haloes 
from the others. It is worth noting that the excess power is 
not simply caused by the presence of these blobs: the wave 
numbers with excess power cover a range from < 100 pc to 
~ 600 pc, and show no sign of convergence on small scales. 
This means there is excess power all the way down to our 
resolution scale of w lOpc, smaller than the size ~ lOOpc of 
the individual over-dense blobs, as well as on scales much 
larger than blobs. Our interpretation is therefore that the 
over-dense blobs are a result of a highly lumpy and turbu- 
lent environment. Below, we will hypothesize that this is a 
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Figure 13. Power spectra of the gas density. For most haloes, the power spectra roughly follow power-law scalings of P p <x k a with a ~ 
on large scales (expected in the limit of negligible pressure and strong shocks) and a <J — 2 (expected in the presence of non-negligible 
pressure and weak shocks) on small scales. Two of the most massive haloes show a clear excess power on small spatial scales (large 
k). These two haloes host over-dense blobs in their outer regions. We interpret this to be a result of a highly lumpy and turbulent 
environment, caused by the biased location and the active merger history of these two haloes. 



result of these haloes biased location (at the knots of the 
cosmic web) and their merging and cooling history (these 
haloes are among those with the largest number of mergers, 
and could cool their gas early on z > 16). 



Figure [L4\ shows the gas velocity power spectra (com- 
pensated by k ). This plot is useful in order to differentiate 
between a Kolmogorov (oc fc _5//3 , or a = 1/3) and a Burgers 
(oc fc -2 , or a = 0) power spectrum. Interestingly, Figure 



14 



shows a flat PS for all 19 haloes, with an average power- 
law exponent of (a) = —0.15 ± 0.019. In other words, the 
gas component for these 19 ACHs is closer to the Burgers 
PS (consistent with earlier findings in Prieto et al.| ( |2011| 
2012)). Figure 15 shows the solenoidal to total (solenoidal 
plus compressional) kinetic energy ratio. In an ideal fully de- 
veloped turbulent fluid, this quantity should take the value 
Rk = 2/3, showing the equal distribution of energy in two 
solenoidal, divergence-free V • v s> k = 0, modes and a sin- 
gle compressional, curl-free V x v c ,h = 0, mode. Figure [15] 
shows that the distribution of energy typically exceeds the 
equipartition value of 2/3; Rk w 0.6 — 0.9. This implies that 
turbulence has not had time to settle into a steady pattern. 



3.5.2 Probability distribution functions of the IHM 

We next examine the probability distribution functions 
(PDFs) of the gas density, temperature, and velocity. 
Throughout this discussion, the PDF is volume-weighted 
(i.e. p(q) gives the probability to find a fluid element with a 
given physical quantity between q and q + dq). 

Figure [16] shows the PDF of the normalized number 
density n/no, where no is the mean number density inside 
the volume analysed (i.e. inside 2_R v ir). Each of the 19 haloes 
shows a distinct peak at ~ 10~ 1 n( ) . Similar peaks can be 



seen in Figure 9 of Prieto et al. (20121; in that paper, the 



peaks show a clear evolution with time, and become more 
pronounced at lower redshifts. As in |Prieto et al.| ( [2012[ ), the 
peaks shown in Figure [16] result from the gas collapsing to 
the central regions of the DM halo. Through this process, 
the gas piles up in a very small percentage of the volume, 
while the low density regions fill the rest of the volume. At 
densities above ~ 10 1 no, the PDF develops a nearly power- 
law shape (the dashed lines in Fig. [16] show a p ower-law 
PDFoc n -3 / 2 for reference). This is explained by 



Kritsuk 



et al. (20111 as a feature from the first stages of a nearby 



isothermal gravitational collapse. 

Two of our 19 haloes developed high-density tails, ex- 
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Figure 14. Gas velocity power spectra (weighted by k 2 ; solid curves). The spectra tend to be very flat. The long-dashed lines show 
power-law fits k 2 P p oc k a in the range 8kpc — 1 < k <80kpc _1 . The average exponent in this range is found to be (a) = —0.15 zt 0.019 
i.e. closer to a Burgers (a = 0) spectrum . 




Figure 15. The ratio between solenoidal and total (solenoidal plus compressional) modes of the kinetic energy. In fully developed 
turbulence, this ratio takes the value 2/3, because the energy is equally distributed among the degrees of freedom (two solenoidal modes 
and one compressional mode). The ratio is found to be 0.6 — 0.9 for our 19 haloes. 
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Figure 16. Number density PDF (solid curves), normalized by the mean density no. The PDFs peak near n ~ 10 _1 no, indicating that 
most of the gas remains at low density, as a small fraction of the gas collapses to high density near the center of the DM halo. Above 
n ru lOnrj, the PDFs develop power— law shapes. The dashed lines show a PDF <x n~ 3 / 2 , characteristic of the first stages of isothermal 
gravitational collapse. Two of the most massive haloes develop a tail at densities above sinl0 4 n/no. These two haloes host massive 
over-dense blobs near their outskirt. 



tending above ~ 10 4 no. These two haloes are the ones that 



exhibit large density fluctuations at small scales (Figure 13 1, 
and develop the compact over-dense blobs (COBs). Such 
high-density tails result from self gravity acting on the gas, 
and have been seen in different contexts in simulations that 
include self-gravity, such as isothermal supersonic turbu- 
lenc e (e.g|Kritsuk et al.|2011[ ), isothermal MHD turbulence 
(e.g. Collins et al. 2012 1 and in a cosmological simulations 



of primordial gas collapsing in a lower-mass ~ 3 x 1O 7 M0 
DM halo ( jPrieto et al.|20T2| ). An interesting feature of these 
two haloes is that they have the broadest number density 
PDFs. This signature is consistent with these two being the 
most shocked and turbulent haloes of the sample. 

Figure |17| shows the gas temperature PDFs for each of 
the 19 haloes. The vertical dashed lines correspond to the 
virial temperatures computed assuming a mean molecular 
weight /i = 1.22 (appropriate for our nearly neutral H+He 
gas). The shape of the temperature PDF can be understood 
as a mixture of cold and hot gas interacting through the 
collapse process. All PDFs show a sharp cut-off above ~ 
10 4 K as a consequence of atomic cooling, and (except for 



the two smallest haloes #1 and #2) the PDFs do not extend 
to temperatures as high as T V - 1T . The PDFs all show cold 
tails, corresponding to uncompressed gas from the IGM. The 
PDFs for haloes #16 and #18 extend close to ~ 10 K; as 
mentioned above, this is due to the activation of a numerical 
temperature floor in the COBs in these haloes. 

Figure [18] shows the PDF of the gas velocity modulus, 
defined as 

V = yj (v x — V x ,cm) 2 + (Vy — Vy^CAl) 2 + (v z — Vz : cm) 2 (19) 

where Vi t CM is the center of mass velocity of the gas in each 
Cartesian direction. For reference, the short-dashed curves 
show Maxwell-Boltzmann distributions assuming a gas tem- 
perature T — Tvir. The velocity PDF resembles features of a 
Maxwell distribution: for most of the haloes it has a power- 
law tail at low velocities that is well appr oximated b y the 
slope PDFoc v 3 , shows a knee at v « -y^/csTvir/m, and 
has an exponential cutoff at higher velocities. For reference, 
the long-dashed vertical lines correspond to the halo circular 
velocity v circ = (GAfvir/^vir) 1 ^ 2 - Interestingly, the peak of 
the velocity PDFs are at high velocities, only slightly below 
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Figure 17. Gas temperature PDF (solid curves). The vertical lines show the virial temperature of each halo. The temperature PDF 
shows a mixture of hot and cold gas, including (i) cool, low-density gas far away from the halo center, (ii) hot, low— density gas, shocked 
during its collapse, and (iii) hot, high— density gas in the central region of the halo. The temperature PDFs are truncated above ~ lO^K 
by atomic cooling. 



(and in some cases coincident with) iWc- This is in con- 
trast with the temperature PDF, which shows that the gas 
cools well below the virial temperature. This indicates that 
a large fraction of the gas falls into the haloes supersoni- 
cally, at speeds near the halo circular velocity. Typical Mach 
numbers of the gas can thus be roughly estimated simply as 
the ratio of the halo circular velocity and the gas sound 
speed A4 = v C irc/c s . Using the spherical collapse model to 
compute v C irc for an ACH of mass M and redshift z, and 
adopting c s = c s (T = 10 4 K), we have 

showing that our haloes should be able to develop a tur- 
bulent environment with Mach numbers A4 ~ 2 — 4, and 
with more massive haloes having higher Mach numbers. We 
computed the Mach numbers inside the virial radius for our 
haloes, and for the most interesting cases, i.e. the RSCs 
and the COBs haloes, we found a Mach number M <4 and 
M > 5, respectively. 



3.5.3 Three-dimensional morphologies 

Figures [T9| and [20| show logarithmic projections of the mass- 
weighted gas density and temperature, respectively, for each 
haloQ There are significant variations in the details of the 
morphologies of our haloes, but all of them are characterized 
by a central overdensity, fed through gas filaments, tracing 
the underlying DM structure. 

The filamentary accretion produces a shocked and tur- 
bulent environment. As mentioned above, the spherical col- 
lapse model predicts that our haloes have Mach numbers 
M ~ 2 — 4 and a direct computation inside the virial ra- 
dius gives a Mach number range of 3 < M < 6 allowing the 
development of supersonic turbulence. One role of the gas 
filaments in this context is to feed the chaotic turbulent en- 
vironment at the central regions of the haloes. As shown in 
Figure |20| (and supported by the average radial profiles in 
Figure |9j) due to the atomic cooling, the overall gas collapse, 
and the accretion through the filaments, both proceed su- 

1 The fourth column in the second row shows two haloes, the 
halo we analyzed (#7) is the one close to the center of the image. 



The formation of disks and SMBHs 17 




-0.5 0.5 1 1.5 2-0.5 0.5 1 1.5 2-0.5 0.5 1 1.5 2-0.5 0.5 1 1.5 2 



log(v/km/s) 



Figure 18. Gas velocity PDF (solid line). In most haloes, a peak is located near the halo's circular velocity (the latter shown by 
long-dashed vertical lines). The short-dashed curves show, for reference, Maxwellian velocity distributions for a temperature T = T v ; r . 



personically. As mentioned in previous sections, two of the 
most massive haloes develop M ~ 1O 7 M0 COBs, which can 
be seen as high density peaks in the bottom row in Figure fl9] 
(#16 and #18). These two haloes present a merger history 
characterized by almost continuous episodes of minor merg- 
ers as can be seen from Table [2] The multitude of mergers 
likely helped to create and maintain a lumpy, turbulent en- 
vironment. 

3.5.4 Turbulence 

The gas in the ACHs have Mach numbers Mi > 1 and they 
are thus susceptible to developing shocks and turbulence 
(see also |Scannapieco fc Briiggcn (2010) for studies of tur- 
bulence in galaxy formation). Our results indeed indicate 
the presence of supersonic turbulence. There are two mech- 
anisms that could naturally drive and maintain this turbu- 
lence: (i)the filamentary accretion of gas in relatively nar- 
row, collimated streams - the infall speed in these streams 
is supersonic, and such streams can collide and create large 
Mach-number shocks, and stir the gas near the core of the 
halo, (ii) mergers with other DM haloes - mergers can like- 
wise create violent shock waves, and stir the gas. 

To elucidate the role of these processes, we first examine 



the polar patterns of mass accretion rate in our haloes. In 
Figure |21| we show the mass infall rate at each point on the 
sphere of radius 0.5i? v i r , placed around the center-of-mass 
of each halo. These polar plots show the mass flux (in units 
of Mq /yr/deg 2 ). Note that negative values correspond to 
inward mass flux (i.e. accretion), and positive values corre- 
spond to outward-directed flows. 

This figure shows sharp peaks, indicating that most of 
the gas inflow into the haloes occurs along relatively nar- 
row accretion streams. We have confirmed that this is true 
quantitatively, i.e. most of the mass inflow occurs in a solid 
angle <C 4-7T. Some haloes (e.g. #5) have multiple streams, 
and some (e.g. #10) have strong collimated outflows as well. 
Most importantly, this figure shows that the two haloes 
in which COBs develop (#16 and #18) have the highest 
mass fluxes by far. Additionally, the streams feeding gas 
into these haloes are the most narrowly collimated (occu- 
pying the smallest solid angle). By comparison, the haloes 
having rotationally supported cores present the lowest mass 
inflows, and the inflow occurs spread out over a much larger 
solid angle (roughly confined, however, to lie within the same 
plane). 

The difference in mass flux among these two populations 
is a factor of ~ 10 1 — 10 2 . These results suggest a causal 
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Figure 19. The logarithm of the projected mass-weighted density. Despite significant variations in their details, all haloes are charac- 
terised by a dense core, and filamentary accretion. Two of the most massive haloes (#16 and #18) develop compact over-dense blobs, 
here seen to be located near the haloes' outskirts. 



connection between the strength of the filamentary accretion 
and the nature of the object formed at the center of the halo. 
Namely, that stronger streams lead to stronger shocks and a 
more turbulent environment, which can produce COBs. By 
contrast, the less collimated, and weaker accretion proceeds 
in a more orderly fashion, which allows the specific angular 
momentum of the gas to be conserved, and a rotationally 
supported core (RSC) to develop. 

Next, we examine the merger histories of the haloes. In 
Table [2] for each halo, we list the total number of merg- 



ers it has experienced, the number of minor and major 
mergers (defined as having mass ratios below and above 
M1/M2 = 1/3, respectively), the redshift and mass ratio 
of each merger, and the first (highest) redshift at which the 
gas in the halo could cool. For our available hydro outputs, 
Zcooi was computed as the z when the gas temperature inside 
7? v ir reaches the T = 10 4 K. Because we do not have all the 
hydro outputs at high z for every simulation, for the z* cases 
in the last column of Table [2] we used the DM-only Simula- 
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Figure 20. The logarithm of the projection of the mass-weighted temperature. These projected images show filaments converging 
towards the central regions of the haloes, producing a shocked and turbulent environment as shown by the values of the structure 
functions (see text) 



tions to get this redshift. In this case z coo i was computed as 
the z when the halo reaches the mass M = M(T vir = 10 4 K). 

This table shows that the haloes with the COBs are 
among those with the highest number of mergers. This is 
mostly due to their being among the most massive haloes. 
However, note that haloes #17 and #19 have comparable 
masses, but experience fewer (i.e. two and zero) mergers, and 
these halo do not develop COBs. This suggests that merg- 
ers are responsible for creating the high mass inflow rates, 
and causing a highly turbulent and lumpy environment that 



leads to strong shocks, and the formation and spatial dis- 
placement of the COBs. Note that the COB haloes have not 
only suffered more mergers, but their gas has also cooled ear- 
lier (see also Fig. [23] below) . This can explain why in these 
haloes the gas is more lumpy (as seen in their density PDFs 



in Fig. 16 1 and can transfer its angular momentum and fall 
in rapidly. Catastrophic transfer of angular momentum and 
rapid gas infall of highly clumped gas to small radii is well 
known to be associated with early cooling in the low-redshift 
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Figure 21. Mass flux across the surface of a sphere with radius 0.5R v i r , centered around each halo. The mass inflows are quoted in 
units of M(j/yr/deg 2 . The haloes that develop compact over-dense blobs (#16 and #18) show the highest inward mass accretion rate 
(several X 10 — 4 MQ/yr/deg 2 ), ~ 10 1 — 10 2 times higher than the accretion rate in haloes with rotationally supported cores (#1, #7 and 
#12). In addition, the mass inflow in the haloes with COBs occurs along the most narrowly collimated streams. This suggests a causal 
connection between the strength of the filamentary accretion and the nature of the object formed at the center of the halo. 



galaxy formation literature (e.g. Weil et aLj 1998 and refer- 
ences therein). 



angular momentum can be conserved, leading to rotationally 
supported cores. 



The same effect can explain the rapid inflow rates we 
find in our two COB-forming haloes. By comparison, in 
haloes that undergo fewer or no mergers, the accretion rate 
remains much lower, spread out over a larger solid angle, and 
more coherent. In this case, a large fraction of the specific 
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Table 2. Merging and cooling histories of the 19 haloes in our sample. The columns labeled Z{ show the rcdshift of each individual 
merger, with the halo mass ratio given in parenthesis. The last column indicates the earliest redshift in which gas cooling in each halo 
was activated. 



Halo Halo mass #Total #Minor #Major zi Z2 Z3 Z4 Z5 zg Z7 z C ool 

(number) Mq mergers mergers mergers (merger 

ratio) 



1 5.35X10 7 - - - _______ 12.36 

2 1.02x10 s - - - _______ 13.94 
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4 DISCUSSION 

4.1 No disks in atomic cooling haloes 

By definition, the ACHs have a virial temperature 
T v i r > W 4 K. These haloes lose part of the thermal energy 
generated through the gravitational collapse process by H 
(and He near ~ 10 5 K) emission line cooling, and reach a 
temperature floor of T ~ W 4 K. This fact can be seen in 
the temperature profiles in Figure |9j which show a constant 
T ~ 10 4 K inside ~ -R v ir- The gas density profiles approxi- 
mately follow a power-law p oc r~ 2 , in-between the profiles 
expected from pressureless cosmological infall p oc r~ 1 5 and 
infall of non-radiating gas through shocks p oc r~ 2 ' 25 ( Gunn 
fc Gott|1972| |Bertschinger|1985 



One of the main motivation for this work was indeed 
to look for large (~ 0.1i? v ir), rotationally supported gaseous 
disks inside the ACHs. This is expected in the simplest toy 



ficient cooling allows the gas to collapse preferentially in the 
direction parallel to its angular momentum vector, while ac- 
quiring increasing rotational support in the plane perpendic- 
ular to this direction - eventually settling to a rotationally 
supported disk. Interestingly, we find that collapse indeed 
occurs preferentially along specific directions, and the gas 
has a less spherical shape than the host DM haloes, haloes 
with masses above ~ 6 x 10 s Mq, 100% present a flattened 
gas distribution, with c/a < 0.5. This is in agreement with 



the earlier results by Romano-Diaz et al. (20111. However, 



models (e.g Mo et al. 1998 Oh fc Haiman 2002 1 , in which ef- 



we find that the minor axis of the gas distribution is not 
aligned with its angular momentum vector - the offset is 
between 30-90 degrees. This invalidates the above simple 
picture. 

Furthermore, we find that none of these flattened gas 
distributions are rotationally supported; their spin parame- 
ters are well below unity. Interestingly, in 3 of our haloes the 
gas in the core does reach Ag as ~ 1, but these rotationally 
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supported cores are found in the lowest-mass haloes, and 
are closer to spherical blobs than to geometrical disks. In 
summary, we have proposed the following requirements for 
the gas distribution to qualify as a conventional "disk" : 

• a low c/a <J0.5 and a high b/a w 1, 

• a high gas spin parameter Ag as w 1 and 

• an alignment (to within ~ 20°) of the angular momen- 
tum vector of the gas with the minor axis of the associated 
ellipsoid. 

We find no disks in our sample of 19 atomic cooling haloes 
that satisfy this set of intuitive criteria. 



4.2 Lack of correlation between gas and DM spin 

Another interesting result from our study is that the direc- 
tions of the angular momentum vectors of the gas and DM 
components are significantly misaligned. The misalignment 
angle 6 > 20° for most of the haloes; in fact, for a third of 
the haloes in our sample, we find 9 > 90°, i.e. the gas and 
the DM components counter-rotate. 

Following the simple scenario in which the baryonic 
matter acquires an average SAM J g as/M gas oc Jdm/Mdm, 
as assumed in |Mo et ah] ( |1998| ) and |Oh fc Haiman| ( [2002] ), 
we should expect a negligible misalignment angle. Three of 
our haloes present alignment to within 6 « 20°. As men- 
tioned above, these haloes have a high spin parameter and 
a high fb, suggesting that they are in a more advanced state 
of collapse compared with the rest of the haloes. One reason 
for the high misalignment angles in most of the haloes could 
be that the gas in these haloes have not yet fully collapsed, 
and the DM and gas components have not had enough time 
to interact and reach a stable, aligned state. As mentioned 



above, Kimm et al. (2011 1 have recently studied the relative 
angular momenta of gas and DM in haloes at lower red- 
shifts, and also found significant differences, due to the very 
different dynamics of the gas and the DM once they cross 
the central region of the halo. It has also been pointed out 
by Angles- Alcazar & Dave (in preparation). 



4.3 Rotationally supported cores vs. compact 
over-dense blobs 

Our simulations show evidence for the formation of RSCs 
and COBs. The haloes in which these two different gas con- 
figurations are found have very different physical features, 
in terms of turbulence, merger history and larger-scale envi- 
ronment. These difference can help us understand the reason 
for having such different primordial objects. 

The radial profiles of the specific angular momenta and 
the spin parameters show that in the 3 haloes where the 
gas component reaches a rotational support at small ra- 
dius, A' w 1 and a L/Lke P ~ 1 inside r w 0.1-Rvir- On 
the other hand, the haloes forming COBs have low gas spin 
parameters, and much flatter spin parameter profiles, indi- 
cating that the gas could efficiently transport its angular 
momentum. As described above, the two types of haloes 
have very different merger histories: Table [2] shows that the 
COB haloes have suffered many more mergers than the RSC 
haloes and that they have cooled their gas earlier (Fig. |23[). 
The early cooling and the multitude of mergers can explain 
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Figure 22. The ratio of the power-law exponents of the longitudi- 
nal structure functions of 2° and 3° order R = £2 /d ■ This ratio 
is expected to be PS 0.74 for fully developed, highly supersonic 
turbulence, and higher for less-developed, weaker turbulence (see 
text). Halo type (i.e. haloes with COBs or RSCs) correlate with 
both this value and mass of the DM halo, hinting at turbulence as 
a possible explanation of the different morphologies of these ob- 
jects. The error bars correspond to the r.m.s among the different 
haloes for each class of objects. 
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Figure 23. The earliest redshift at which atomic cooling is ac- 
tivated in each of the 19 haloes. The haloes forming COBs are 
among the most massive at 2 = 10 and are the first to cool. 



the fact that the gas in the COB haloes have a much more 
clumpy distribution (Fig. 



161, which may explain the ef- 



ficient transfer of angular momentum, and the rapid and 



narrowly collimated gas inflows (Fig. 21 1 in these haloes - 
creating strong shocks and driving strong turbulence in the 
central regions, eventually forming the COBs. 

In order to distinguish haloes with a diagnostic that is 
explicitly related to turbulence, we have measured the ve- 
locity structure functions of the gas within the inner 0.5i? v ir- 
The structure functions of order p are defined as: 
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M DM =5.354x1 7 M SLln . RSC (#1 ). 



M DM =9.220x1 B M sun . COB (#1 6). 
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Figure 24. Maps of the projected DM surface density, in arbitrary units, to illustrate the large-scales cosmic environments of one of 
the haloes developing a rotationally supported core (RSC, left panel) and a compact over-dense blob (COB, right panel). The RSC halo 
is located at a typical point in a filamentary region, and the COB-forming halo resides at the center of a dense knot in the cosmic web. 
The merging histories of these two haloes are accordingly also very different (see Table [2| . 



S P {1) = (\u(x + 1)- u{x)\ p ) oc Z c(p) 



(21) 



where the velocity component it is parallel (longitudinal 
structure function) or perpendicular (transversal structure 
function) to the vector I, I =j I j, and the spatial average 
is over all values of the position x (for more details, see 
Padoan et al. (20041). The ratio of the power-law exponents 
of the longitudinal structure function of 2° and 3° order, 
R — C2/C3! has been shown to depend on the strength of 



turbulence (Padoan et al. 20041. In particular, for highly 
developed supersonic turbulence, |Padoan et al.| ( |2004| finds 
the general scaling 



C(2) 
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+ 1 - (1/3) 



(P/3) 



(22) 



which for p = 2 yields a value of 0.74. We have computed this 
ratio for each of our 19 halos, and show the average values, 
with r.m.s. errors, for the three types of halos separately in 
Figure [22] Interestingly, the two types of halos with RSCs or 
COBs have significantly different ratios, -Rrsc = 0.90 ±0.04 
versus i?coB = 0.79 ± 0.01, which can be attributed to 
the different turbulent conditions. The value for the COB- 
forming halos is close to the Boldyrev-type (i.e. highly super- 



sonic, fully developed) turbulence. From Figure 2 in Padoan 



et al. (20041, we see that the higher value of R — 0.9 for the 



halos with RSCs corresponds to less supersonic, and less 
well-developed turbulence (but still above the value for Kol- 
mogorov subsonic turbulence). Perhaps most interesting is 
the clear correlation between the halo mass and structure 
function ratio, as can be seen in Fig. |22| Without further 



simulations and analysis, it is unclear what causes this cor- 
relation, but we speculate, consistent with our interpreta- 
tions above, that it is due to the different merger histories, 
which in turn set the different turbulent conditions, which 
in turn create different types of objects. We will explore this 
hypothesis in more detail in a forthcoming publication. 

As mentioned above, the mergers suffered by the pro- 
genitor halo is not the only possible source of turbulence. 
A second possible source is rapid gas accretion through 
filaments, following the underlying DM structures. These 
two mechanisms can, in principle, both simultaneously con- 
tribute to driving turbulence. In particular, we found that 
the gas accretion rate for the COB haloes is ~ 10 — 10 2 times 
higher than the one presented in the RSC haloes, suggesting 
that this second mechanism must play an important role. 

Finally, it is interesting to ask whether the larger-scale 
environments of the haloes may affect their turbulent prop- 
erties, and the type of object that forms the haloes. Figure 
[24] shows maps of the projected DM surface density around 
the location of one of the RSC (left panel) and COB (right 
panel) haloes. The differences are evident: while the RSC 
halo is located at a typical point in a filamentary region, the 
COB-forming halo resides at the center of a dense knot in 
the cosmic web. These DM density maps support the idea 
that the more numerous mergers and more rapid accretion 
flows found in the COB haloes are connected to its larger- 
scale cosmic environment. 
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5 SUMMARY AND CONCLUSIONS 

In this study, we have presented the largest statistical cos- 
mological hydro-simulation study of atomic-cooling haloes 
(ACHs). We studied the physical properties of 19 isolated 
haloes at z = 10, covering a total comoving volume of ~ 2000 
Mpc 3 . The simulations reach a spatial resolution of Ax ~ 10 
(proper) pc. 

We presented a comprehensive analysis of the gas and 
DM physical properties, including baryonic mass fractions, 
gas and DM shapes, gas and DM spin parameters, as well as 
the radial profiles, power spectra, and probability distribu- 
tion functions of various halo and gas properties. Our main 
goal was to study and understand the formation of any large 
galactic disks among ACHs, which may be stable and avoid 
fragmentation, suitable for forming SMBHs. In particular, 
our aim was to assess whether disk formation was common 
among these haloes. 

Our main conclusions can be summarised as follows: 

(i) None of our 19 haloes has formed a rotationally sup- 
ported geometrical disk. While haloes with masses above 
~ 10 9 Mq have flattened geometrical shapes, they do not 
have significant rotational support. Conversely, we find three 
lower-mass haloes in which a rotationally supported core de- 
velops, but the shape of this core is closer to spherical than 
to a disk. Nevertheless, these cores are candidate sites for 
SMBH formation. 

(ii) The gas spin parameter Ag as is higher than the DM 
spin parameter Ad M for some of our haloes above ~ 1O 8 M . 
This result is in superficial agreement with the simple pic- 
ture of galactic disk formation, in which specific angular mo- 
mentum is conserved, and the gas spin parameter increases, 
as the gas collapses inside the DM halo. However, we find 
that the magnitude of the gas and DM spins have little corre- 
lation, and their directions are also significantly misaligned. 
This indicates that the gas angular momentum does not fol- 
low the above simple picture. 

(iii) Two of our most massive haloes developed very 
massive compact over-dense blobs, with masses of 
~fewxlO 7 M . The COBs form near the center of the haloes 
by strong shocks at the intersection of converging filamen- 
tary flows. However, unexpectedly, the blobs have large 
residual velocities, and they migrate to the haloes' out- 
skirts. The blobs are self-gravitating, and accrete rapidly 
(~ 0.5M Q yr _1 , measured at a distance of r « 100 pc from 
their center). These COBs show no sign of fragmentation at 
our resolution, and are alternative candidate sites for SMBH 
formation by direct collapse. 

(iv) The intra-halo medium is turbulent and supersonic 
(with Mach numbers M = 2 — 4). To demonstrate this we 
have computed structure functions and showed that for all 
haloes they are all consistent with structure functions above 
the Kolmogorov regime. 

(v) The COB formation is associated with a highly tur- 
bulent, supersonic (Mach numbers M > 5) and lumpy en- 
vironment due to a high mass accretion rate and an active 
merger history. This is likely related to their location at 
dense knots of the cosmic web. On the other hand, the rota- 
tionally supported cores are associated with haloes that are 
less turbulent, undergo mergers much less frequently, and 
have a lower gas accretion rate. They are located midway 
along filaments in the cosmic web. 



We conclude that the angular momentum transport and 
the fate of gas in z = 10 ACHs is intimately related to the 
turbulent conditions, which in turn is determined by the 
merging history and larger-scale environment of the halo. 
We hypothesize that haloes in highly biased regions have 
more numerous and earlier mergers, more lumpy and tur- 
bulent gas distributions, and more efficient transfer of an- 
gular momentum. The gas in such haloes can be accreted 
rapidly, along narrow streams, into a very dense central ob- 
ject, whose formation involves strong shocks. On the other 
hand, haloes in less biased regions have fewer mergers, and 
a less turbulent intra-halo medium, where gas inflow is more 
coherent and weaker overall. In such haloes, a larger fraction 
of the angular momentum can be preserved, and thus even- 
tually a rotationally supported configuration can be reached. 

Our results are based on simplified simulations. Since we 
assume the gas has primordial composition and H2-cooling 
and prior star-formation in the haloes have been suppressed, 
our simulations can be regarded as numerical experiments. 
They may nevertheless be an accurate description of a small 
subset of atomic cooling haloes, exposed to strong enough 
Lyman- Werner radiation to prevent prior star-formation. If 
verified in a larger sample of haloes and with additional 
gas physics to account for metals, star-formation and feed- 
back, our results will have implications for observations of 
the highest-redshift galaxies and quasars with JWST. 
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